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ABSTRACT 

We calculate the structure of accretion disc around a rapidly rotating black hole 
with a super-Eddington accretion rate. The luminosity and height of the disc are re- 
duced by the advection effect. In the case of a large viscosity parameter, a > 0.03, the 
accretion flow strongly deviates from thermodynamic equilibrium and overheats in the 
central region. With increasing accretion rate, the flow temperature steeply increases, 
reaches a maximum, and then falls off. The maximum is achieved in the advection 
dominated regime of accretion. The maximum temperature in the disc around a mas- 
sive black hole, M = IO^Mq, with a = 0.3 is of order 3 x 10^ K. Discs with large 
accretion rates can emit X-rays in quasars as well as in galactic black hole candidates. 

Key words: accretion, accretion discs - black hole physics - hydrodynamics - rela- 
tivity - radiation mechanisms 



1 INTRODUCTION 

Black hole accretion discs are believed to be the mechanism 
of energy release in active galactic nuclei and in some X-ray 
binaries. Their luminosity is constrained by the Eddington 
limit, Le = 4:TTGMmpc/aT {M is the black hole mass), and 
the accretion rate, M, is usually expressed in the Eddington 
units, m — Mc^/Le- The luminosity becomes equal to the 
Eddington limit at rhcr ~ , where r) is the radiative 
efficiency of the disc. Accretion rates of order likely 
to feed the huge energy release in quasars, and m ~ rhcr is 
implied in many models of the soft X-ray excess observed in 
quasar spectra (e.g. Czerny & Elvis 1987; Dorrer et al. 1996; 
Szuszkiewicz, Malkan & Abramowicz 1996). 

According to the standard model (Shakura 1972; 
Shakura & Sunyaev 1973), the effective temperature of the 
disc surface scales as Te// oc (m/Af)^''^. A special feature 
of discs with m ~ rhcr is that their temperature can be 
much greater than T^jj provided the viscosity parameter 
Q is large enough (Shakura & Sunyaev 1973). In the inner 
region of such a disc, the inflow time-scale is shorter than 
the time-scale for relaxation to thermodynamic equilibrium 
and the accretion flow is overheated. The overheating might 
be strong enough for a transition from the standard radia- 
tion dominated disc to the hot two-temperature ion-pressure 
supported regime of accretion proposed by Shapiro, Light- 
man & Eardley (1976). This transition has been discussed in 
recent works (Liang & Wandel 1991; Artemova et al. 1996; 
Bjornsson et al. 1996). 

The standard model can be a good approximation only 
if m < rhcr- When the accretion rate approaches rhcr the 
main assumption of the model, that the viscously released 



energy is radiated away locally, becomes unadequate. Inside 
some radius rt, the produced radiation is trapped by the flow 
and advected eventually to the black hole instead of being 
radiated away. This radius can be estimated by taking the 
internal energy density, e, and the density, p, of the disc 
from Shakura & Sunyaev (1973). Then one finds that the 
radial flux of internal energy = Me/ p exceeds the total flux 
of radiation emitted between r and 2r at 



2GM/c? is the gravitational radius. For m 



r,^hrr,^l-(^^Y\ (1) 
where 

rhcr = 12 equation (1) yields rt ~ T-lrg. A substantial por- 
tion of the released energy must be swallowed by the black 
hole even when rh < rhcr- On the other hand, the relative 
height of the Shakura-Sunyaev disc, H/r, equals rh/27 at 
the maximum. So, advection becomes essential before the 
accretion flow becomes quasi-spherical, and this efi'ect can 
be investigated in an extended version of the standard model 
retaining the vertically integrated approximation. The cor- 
responding set of equations was proposed by Paczyiiski & 
Bisnovatyi-Kogan (1981). Numerical solution of these equa- 
tions gives adequate description for the inner transonic edge 
of the accretion flow and shows that the disc remains rel- 
atively thin at moderately super-Eddington accretion rates 
(Abramowicz et al. 1988; Chen & Taam 1993). This made 
natural the extension of the standard model to the super- 
Eddington advection dominated regime, called "slim" accre- 
tion disc. 

All the models of the super-Eddington slim disc em- 
ployed the pseudo-Newtonian approximation to the black 
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hole gravitational field (Paczynski & Wiita 1980) which is 
good only in the case of Schwarzschild black hole. Recently, 
the equations of relativistic slim disc has been derived by 
Lasota (1994, see also Abramowicz et al. 1996) and applied 
to another type of advection dominated accretion flow, op- 
tically thin and hot. 

In the present paper, we solve the relativistic equations 
of the super-Eddington disc. In the case of a rapidly ro- 
tating black hole relativistic effects become especially im- 
portant. The sonic radius is typically inside the ergosphcrc, 
close to the black hole horizon. The released power steeply 
increases when the hole spin approaches its extreme value, 
a* = JcjM^G — > 1 ( J is the angular momentum of the black 
hole). In this paper, wc calculate the disc structure for the 
case a* = 0.998 and compare it with the case of non-rotating 
black hole a, = 0. 

We pay particular attention to the case of a large vis- 
cosity parameter, a. In the previously investigated pseudo- 
Newtonian models, the supor-Eddington disc was assumed 
to be in thermodynamic equilibrium which is a good ap- 
proximation in the case of small a. We find that a signifi- 
cant deviation from the equilibrium occurs when a ^ 0.03. 
Then the flow overheats in the central region. We find that 
the temperature of the fiow steeply rises at rn > rhcr and 
reaches a maximum in the advection dominated regime of 
accretion. The maximum temperature is especially high in 
the case of a rapidly rotating black hole. 

In next Section, we write down the equations of the rel- 
ativistic advcctive disc. The equations describe a radiation 
dominated fiow and neglect the thermal pressure of the ac- 
creting plasma. Their solution yields the radiation density, 
plasma density, and height of the disc within the equipar- 
tition radius of the standard model, r^q, at which the ra- 
diation pressure equals the plasma pressure. The transition 
from the standard model to the advective accretion occurs 
at rt req- We describe the numerical method, present 
the solution, and discuss the structure of the disc in Sec- 
tion 3. Then, in Section 4, we calculate the temperature 
of the accreting plasma. Finally, we check that inside req 
our radiation-dominated models never approach the ion- 
pressure supported regime of accretion. The results are sum- 
marised in Section 5. 



2 RELATIVISTIC DISC EQUATIONS 

In the description of the gas flow around a Kerr black hole 
we use the Boyer-Lindquist coordinates = {t, r, d, (p). The 
metric tensor gij of the Kerr space-time is given, e.g., by Mis- 
ner, Thorne & Wheeler (1973). The disc is assumed to lie at 
the equatorial plane of the Kerr geometry. The four- velocity 
of the accreting gas has components — {u* ,u'^) in 
the Boyer-Lindquist coordinates, and = 0. The angular 
velocity of the gas rotation equals Q = uf jx^ and Lorentz 
factor measured in the frame of local observers with zero 
angular momentum equals 7 = «*(— g**)~^/^. 

To obtain a simple model of the relativistic disc with ad- 
vection we employ the vertically integrated (slim) approxi- 
mation. The slim disc equations are written for the following 
thermodynamic quantities; surface rest mass density E, sur- 
face energy density V (which includes both the rest mass 
energy Ec^ and the internal energy 11), and the vertically 



integrated pressure P. The dimcnsionlcss specific enthalpy 
is defined as /i = (?/ -|- P)/Ec^. The equation of state in the 
radiation dominated disc is IT w 3P. F"*" denotes the sur- 
face rate of viscous heating, and F~ denotes the radiation 
fiux radiated from both faces of the disc. All the thermody- 
namic quantities and both fluxes F~ , P"*" are measured in 
the comoving frame. 

The main equations express the conservation laws for 
barion number, energy, angular momentum, and radial mo- 
mentum. We include in the equations the inertial mass as- 
sociated with internal energy accumulated in the flow (Be- 
lohorodov, Abramowicz & Novikov 1997). 



i) Barion conservation 
2'KrcuY. = -M. 

ii) Conservation of angular momentum 



d_ 
dr 



(2) 



(3) 



where v is the kinematic viscosity, and a^^ is the shear, 

r 1 rr / 77 3 

lii ) First law of thermodynamics 

where ^ « 1 is a numerical factor accounting for non- 
homogeneity of the disc in the vertical direction, hereafter 
we set ^ = 1. The surface heating rate is given by 



P+ = 2i.EM cT2c^ = ig" 



iv) Radial Euler equation 



flw [-9 



1 d_ 

2 dr 



dr 



^ffV (n-n+) {n-n-^) 



dP F+i 



c^E/x dr c^E/u 
where are the Keplerian angular velocities. 



(5) 



>(2r/rg)V2±o- 



The set of disc structure equations becomes closed when 
the viscosity, u, and radiative cooling, F~ , are specified. The 
standard a-prescription for viscosity is u = aCeH, where a 
is a constant, Ca = c{P/Uy^^ is the isothermal sound speed. 
The half-thickness of the disc, H, should be estimated from 
the vertical balance condition. 

v) Vertical balance 

Near the black hole, relativistic effects become impor- 
tant and the tidal force contracting the disc in vertical di- 
rection depends on fi. At = fi^, the vertical tidal acceler- 
ation in the comoving tetrade equals (e.g. Riffert & Herold 
1995) 



zrg{r'^ - a,rg^/W~r + 0.75alr 



r^(2r^ — 3rgr + a-^rg^^'lrgr) 



3 Kr 



(6) 



where z = ^/gee (^ — 7r/2) is the height in the disc, J is 
a relativistic correction factor becoming unity aX r ^ rg. 
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Then the typical half-thickness of the disc can be estimated 
as 

More accurate expression accounts for a deviation AQ — 
f2 — of the gas rotation from Keplerian (Abramowicz, 
Lanza & Percival 1997). We will consider only the disc region 
outside the sonic radius where the correction to estimation 
(7) connected with Afi does not exceed several percent, and 
hereafter we use this estimation. 

vi) The radiative losses 

The time-scale for photon diffusion from the interior 
of the disc to its surface equals to = Htq/c where To = 
aT'S/2mp. This is a typical leaking time for the radiation 
trapped inside the disc, and the radiative losses can be writ- 
ten as 



F-=x 



rripC n 



(8) 



where x ~ 1 is a numerical factor. In the standard model 
with the vertical structure governed by the radiation diffu- 
sion this factor equals 2/\/3. In the advection dominated 
region, to exceeds the inflow time-scale ta and a detailed 
treatment of the stationary diffusion is not relevant: ta <tD 
means that the gas accretes faster than a stationary vertical 
distribution of the trapped radiation could establish. How- 
ever, the estimation (8) gives the wright limit F~ <C F^ at 
r <C n. Indeed, when ta >C to we have 11 w F'^ta, hence 
F~ jF^ » ta/tn The radiative losses are "turned off' 
inside rt, and the exact value of F~ is not important for hy- 
drodynarnical behaviour of the flow. A detailed calculation 
oi F~ would require 2D simulation of the radiation diffu- 
sion in the disc with a specified vertical distribution of the 
viscotis energy release. 

For dcfiniteness we hereafter choose x = 2/\/3 in equa- 
tion (8). 

vii) Global energy conservation 

The luminosity of the disc is related to the accretion 
rate by (Beloborodov et al. 1997) 



= -^JutF-rdr^Mc'(l+„s^y 



(9) 



where index "s" refers to the inner transonic edge of the 
disc. The advection effect results in that the luminosity is 
less than the total power released in the disc, 



oo 



utF rdr. 



The power "swallowed" by the black hole equals Ladv = 
L+-L-. 



3 NUMERICAL SOLUTION 

To solve rmriierically the disc structure equations, we choose 
three independent variables f2,Cs, and C = IJ-i^v^ra^^ + 
Muy, /2t:). Prom equations (3), (4) we have 



dr 



P_^_^d_ 
dr dr 



= )^ 



(10) 
(11) 



Expressing u,p,F^,P,Yl in terms of f2,Cs,C we get two 
differential equations for the three unknowns. The third 
differential equation is the radial equation (5). We solve 
the equations (5,10,11) with external boundary conditions 
Cs°'"*,n°'"*,C°"* at a radius r^ut such that n < rout < req. 
The values of c°"* and n°"* = n+ are taken from the rel- 
ativistic version of the standard model (Novikov & Thorne 
1973; Page & Thorne 1974) with the corrected vertical bal- 
ance (Riffert & Herold 1995). The value of C°"* is the eigen 
value of the problem, see below. 

We use the relaxation technique. As first approxima- 
tion we assume f2(r) — n^(r) and solve equations (10,11) 
for Cs-X- Substituting the obtained solution to the radial 
equation (5), we calculate Q'{x) needed to fulfill this equa- 
tion. Then we perform the relaxation step changing f2(r) to 
Q{r)-\-€5fl{r) where 5Q — Q' — Q and e < 1 is the relaxation 
parameter. With new fl{r) we solve again equations (10,11), 
and so on until 5Q converges to zero at all r. This method 
is numerically unstable unless some smoothing procedure is 
used to suppress numerical oscillations in (5f2(r). We have 
found the needed procedure as a combination of two-point 
and three-point smoothing, and determined Af2 = Q — 
with an accuracy of ~ 0.1 %. 

This technique allows to obtain a solution of equations 
(5,10,11) for any given C'"'*- Then we adjust C°"' to fulfill 
the regularity condition at the sonic radius, r^, and find the 
transonic solution. The regularity condition can be written 
as N = D = where A'^ and D are the numerator and 
denumerator in the radial equation with explicitly expressed 
derivative of velocity, du^ /dr — N/D. In this way, we find re 
with an accuracy of ~ 1 % and calculate the disc structure 
at r > Ts (inside the gas is almost in free fall). 

Example solutions outside rs are shown in Pig. 1 for the 
cases of a, = and o* = 0.998 (the represented magnitudes 
are independent of the black hole mass). In these models 
a = 0.1. To a good accuracy, solutions with different a' can 
be obtained from the solutions with a = 0.1 using a simple 
scaling law: E' = Y:a/a' , P' = Pa/a' , with the same H{r), 
AQ{r), and Cs{r). Essential deviation from this law appears 
only when a ^ 1. Then the sonic radius moves outwards 
being beyond the radius of the Keplerian marginally stable 
orbit, rms- At small a, rs < rms- 

The relativistic version of the standard model is shown 
by the dashed lines in Pig. 1. The accretion fiow deviates 
from this model at the trapping radius, rt- Inside rt, the 
condition Me? ^ > 2Tvr'^F~ takes place which means that 
the radial fiux of internal energy exceeds the local radiative 
losses, i.e., the flow is advection dominated. The trapping 
radius is practically independent of a, and the dependence 
on rh is given in Pig. 3. From Fig. 1 one can see that: 1) the 
advection effect reduces the height, H, below the standard 
value, so that the relative disc height, H/r, is kept modest 
even at m ^ lOrhcr, 2) the Thomson optical depth of the 
disc. To, has a minimum of ~ 100 (a/0.1)^^ at rh ~ ihcr, 
and 3) the maximum deviation from Keplerian rotation is 
~ 20 %. 
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Figure 1. The structure of the disc around a black hole with a* = (a-c) and a, = 0.998 (d-f) in the case of viscosity parameter 
a = 0.1. The upper panels show the relative height of the disc, the panels in the middle show the Thomson optical depth of the flow, 
and the bottom panels show the deviation from Keplerian rotation. All these quantities are independent of the black hole mass. Curves 
with different numbers correspond to different a<;cretion rates: 1 - m = rhcr {ihcr = 17.5 in the case of o* = and mcr = 3.11 in the 
case of a* = 0.998), 2 - m = 100, 3 - m = 1000. The dashed lines display the relativistic standard model. 



The critical accretion rate rhcr = 17.5 in the case of 
a* = and rhcr = 3.11 in the case of a, = 0.998. At m > 
rhcr a substantial fraction of the total energy released in the 
disc is advected into the black hole, see Fig. 2. 

Note that in the standard disc with rh < rhcr the den- 
sity, n = 'S/2Hmp, scales as rh~^. In the advection domi- 
nated regime (m > rhcr, r <C rt), the density scales as rh. 



It means that at m ~ rhcr there must be a minimum of the 
density. At this minimum a strong overheating of the flow 
is possible as discussed in next Section. 
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1 10 100 1000 

m=McVLj; 



Figure 2. The fraction of tlie total power released in the disc 
that is eventually advected into the black hole (see item vii of 
Section 2). Curves 1 and 2 correspond to a* = and a* = 0.998 
respectively. 

4 OVERHEATING IN THE CENTRAL 
REGION 

An accretion disc can be considered as a two-fluid flow 
"plasma + radiation". In the case of a small viscosity, 
the disc is dense and the plasma and radiation are close 
to thermodynamic equilibrium at a common temperature 
T ~ Teff = (w/ttr)^''* where w = 3P/2H is the radiation 
density calculated in Section 3, and ar = 7.56 ■ 10^^'' is the 
radiation constant. [Note that Tg// inside the disc differs 
from the surface effective temperature T^ff = (F~/2cr)^/'' 
where cr = arc/4 is the Stefan-Boltzmann constant.] How- 
ever, when a > 0.03 and rh approaches rhcr a strong devi- 
ation from the equilibrium occurs in the inner region of the 
disc. In this case, the accreting plasma has a low density 
and therefore, a low emission capability. As a result, the 
plasma docs not manage to reprocess the released energy 
into Planckian radiation, and the balance between heating 
and hcat-rcrnoval through radiative cooling maintains a tem- 
perature T ^ Tf.ff. Then the radiation is concentrated in 
the Wien peak of temperature T and its density, w, is much 
loss than the corresponding Planckian value, Wpi = arT*, 
i.e., the flow is far from the black body state. 

a) Radius of the overheated region 

In the outer "black body" region of the disc, the free- 
free emission capability, w//, of the plasma with T = Teff 
exceeds the heating rate, w'^ = /2H. The boundary r* 
between the black body and the overheated regions can be 
estimated from the condition 




1 10 100 1000 

m=McVL(. 



Figure 3. The radius of the overheated region r* versus m in the 
case of « = 0.3, M = 10* M©, a* = 0.998. For comparison the 
trapping radius, rt, is also plotted. 

w+=Wff{Teff) at r = r*, (12) 
Wff = 1.6- W'^'^n^Vr. 

At exact thermodynamic equilibrium, free-free emission 
would be balanced by free- free absorption, and Wff can be 
written as wff = cnaffWpi. off introduced in this way has 
the meaning of an effective cross section for absorption of the 
Planckian radiation; it is w 4 times the Rosseland averaged 
Off- 

Radius r* is usually estimated from the condition that 
the effective optical depth of the flow r, = (tqt//)^/'^ — 1 
where Tff = Hnaff (Rosseland (J// is also often taken in- 
stead of iJff). For the standard disc this condition is equiv- 
alent to equation (12). In that case, the radiation density 
inside the disc equals w « w'^tD where to = tqH/c is the 
diffusion time-scale. Then equality (12) with w ~ Wpi gives 
r* « 1, and the decoupling of w below Wpi happens if r, < 1. 
The condition t* < 1 means that a deviation from thermo- 
dynamic equilibrium in the standard disc occurs if the bulk 
of radiation diffuses out without absorption and re-emission. 

This condition, however, is not relevant in the advec- 
tion dominated regime. As the radiation is advected rather 
than escape, the time-scale for free-free absorption should 
be compared with the inflow time-scale rather than with 
the diffusion time-scale. Then the condition r, < 1 should 
bo replaced by r* < {to/taY^^ ~ rt/r that follows from 
equation (12). 

We have calculated r* from equation (12) for the case of 
a = 0.3, M = 10* Mq, a, = 0.998. In Fig. 3, r, is shown ver- 
sus rh and compared with the trapping radius, rt- One can 
see that the entire advection dominated region is overheated 
at rh < 300. 
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h) Cooling rate inside r* 

The overheated plasma inside r* is cooled by Comp- 
tonized free-free emission with rate 



Wcool 



Awff, 



(13) 



where A{n, T) is an amplification factor due to Compton 
upscattering of photons emitted at low energies hv < kT. 
This factor is given by (e.g. Rybicki & Lightman 1979) 



3 2 

^ = 1 + - In Xcoh, 



(14) 



where Xcoh = hi^coh/kT is the dimensionless photon en- 
ergy below which the absorption time-scale is shorter than 
the time-scale for shifting in frequency due to Comptoniza- 
tion. Xcoh is determined by the equation aff{xcoh) = SOar, 
where 9 = kT/m^c^, and (7ff(x) = 1.8-10"^^ In(2.25/a;)(l - 
Expression (14) for the Compton en- 
hancement factor assumes that photons with x > Xcoh up- 
scatter to the Wien peak before they can escape the disc 
or advect to the black hole. The time-scale for upscatter- 
ing to the Wien peak is tc ~ ln{x~^fj /Sdnarc. Hence, 
tc/to ~ ^^i^coh) /y where y = 4:6tq is the Kompaneets' 
parameter. In the considered situation, y ^ 1 and photons 
with X > Xcoh do comptonize to the Wien peak before they 
can escape. Inside the trapping radius, tc should be com- 
pared with ta- We check in the calculated models that tc is 
shorter than ta as well. 

c) H eating = cooling balance 

Everywhere in the disc the plasma temperature is deter- 
mined by the heating=cooling balance w"*" = Wcooi- Inside 
Tt, Wcool is given by equation (13). Outside r,, the cooling 
rate is equal to a difference between free-free emission and 
absorption, being proportional to a small deviation of Wpi 
from w. Wcool = Wff{l — w/wpi). The transition at r ~ r* 
can be smoothly described by the following interpolation for 
the heating=cooling balance. 



+A[l 

Wpl \ Wpl 



w 

Wpl 



(15) 



In the limit w ~ Wpi this equation yields w"*" — ^//(l — 
w/Wpi) while in the limit w Wpi it transforms into = 
Awff. 

d) Results 

Using the balance (15), we have calculated the temper- 
ature in the disc around a massive black hole M — W^Mq 
with spin a* = and a, — 0.998. The results are shown 
in Fig. 4,5 for the case of a = 0.3. The strong overheating, 
T ^ Tcf f , occurs inside the radius r, . The maximum tem- 
perature is reached in the innermost region. It is ~ 3 • 10* 
K in the case of a rapidly rotating black hole and ~ 10^ 
K in the case of a Schwarzschild black hole. In Fig. 6,7 we 
show the dependence of the maximum temperature on the 
accretion rate in the cases of a = 0.03 and a = 0.3. 

A question of interest is whether the protons can be 
thermally decoupled from the electrons at such tempera- 
tures. The time-scale for Coulomb energy exchange between 
protons of temperature Tp and electrons of temperature Tc 




Figure 4. The temperature of the disc with a = 0.3 in the case 
of a, = 0: 1 - m = 10, 2 - m = 100, 3 - m = 1000. The dotted 
Hncs show the effective temperature, Tgff, of the radiation field 
in the disc. 




100 
rcVGM 



Figure 5. The temperature of the disc with a = 0.3 in the case of 
a, = 0.998: l-m = l,2-m=10, 3-m = 1000. The dotted lines 

display the effective temperature, Tcff, of the radiation in the 
disc. The heavy line shows the temperature for which = Ucoh 
at B = Bcq (see the text). 
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1000 



Figure 6. The maximum temperature in the disc as a function 
of the accretion rate for the case of a* = 0: 1 - a = 0.03, 2 - 
a = 0.3. Tgf f at the maximum is displayed by the dotted lines. 




10 

m=McVLj, 



Figure 7. The maximum temperature in the disc as a function 

of the accretion rate for the case of a* = 0.998: 1 - a = 0.03, 
2 - a = 0.3. Tgff at the maximum is displayed by the dotted 
lines. Above the dashed line, ub exceeds I'coh in the inner hottest 
region of the disc. 



is given by (Landau & Lifshitz 1981) 

fcTe 1 



TT JTlp / kTe \ 

2 rrie KnieC^ J 



n3/2 



In A tkjtc 



12.5- 



whcrc In A 



20 is the Coulomb logarithm. Even in the 



most hot models with a = 0.3 and a* = 0.998, tep <C ta 
and the accreting plasma can be well described in the one- 
temperature approximation Te w Tp w T. 

Then we have compared the thermal plasma pressure, 
Pgas = 2nkT, with the radiation pressure, prad ~ w/3. Even 
in the hottest region of the disc with a = 0.3, the plasma 
pressure does not exceed ~ W~^Prad- So, the calculated ra- 
diation dominated models are self-consistent. Note, that ac- 
cretion flows with a ^ 1 are very unlikely (it would imply a 
turbulence of scale H with sound speed). Anyway, such flows 
could not be much hotter because at higher temperatures 
additional cooling mechanisms become important: Comp- 
tonization of cyclotron radiation and collective waves in the 
plasma, as noted by Shakura & Sunyacv (1973). Cyclotron 
radiation and collective waves are additional sources of soft 
photons which can be upscattered to the Wien peak and cool 
the plasma efficiently. Roughly, the cyclotron source starts 
to be important when the gyrofrequency, i^b = eB /2nrncC, 
exceeds Vcoh (a more detailed treatment including higher 
cyclotron harmonics is given in Gnedin & Sunyaev 1973). 
In a similar way, Comptonization of collective waves can be 
important when Upi ^ Vcoh-, where Vpi = [ne^ /■KmeY^'^ is 
the plasma frequency. In case the magnetic field is compara- 



ble to the equipaxtition value 



^l/2 



the gyrofre- 



quency exceeds Vcoh in the disc sooner than the plasma fre- 
quency does. In Fig. 5, the heavy line shows the temperature 

at which the gyrofrequency equals Vcoh- The corresponding 
maximum temperature in the disc is shown by the dashed 
line in Fig. 7. 



5 CONCLUSIONS 

We have investigated the relativistic accretion disc around 
a Kerr black hole with a super-Eddington accretion rate. 
The disc has a modest relative height, H/r < 0.4, up to 
r'n Wrhcr, and a luminosity near the Eddington limit. The 
bulk of the energy released in the inner region of the disc is 
advected into the black hole. The angular velocity of rota- 
tion deviates from Keplerian up to 20 %. We paid particular 
attention to the case of a large viscosity parameter a > 0.03. 
In this case, the accretion flow deviates from thermodynamic 
equilibrium and overheats in the central region. The hottest 
flow has accretion rate rh ~ Srhcr- We have calculated the 
maximum temperature in the disc around a massive black 
hole, M = 10* Mq, with a = 0.3. For a Schwarzschild black 
hole it equals ~ 10^ K. For a rapidly rotating black hole, the 
maximum temperature is of order ~ 3 • 10® K. It far exceeds 
the effective temperature corresponding to thermodynamic 
equilibrium. However, it is not large enough for thermal de- 
coupling of the protons and transition to the ion pressure 
dominated regime of accretion. 

The employed vertically integrated model gives approxi- 
mate characteristics of the bulk of the gas flowing in the disc, 
and docs not describe physical conditions in the upper layers 
(in the "skin" of the disc) where the spectrum of emerging 
radiation is formed. To evaluate the spectrum of the super- 
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Eddington disc a more detailed treatment of the vertical 
energy transfer is needed which is likely to demand a 2D 
simulation of the advective region. Besides, the disc spec- 
trum may be strongly affected by a corona activity above 
the surface. It is clear, however, that the surface of the flow 
with large a and a largo accretion rate must be much hotter 
than the effective surface temperature, T^/f = {F~ /2ay^*, 
and an activity of the corona may make the spectrum only 
harder. X-ray emission of the innermost region of the super- 
Eddington disc may help to reconcile the observed quasar 
spectra with accretion disc models. 
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